Dynamic Jahn- Teller effect in electron transport through single Ceo molecules 
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Scanning tunneling spectra on single Ceo molecules that are sufficiently decoupled from the sub- 
strate exhibit a characteristic fine structure, which is explained as due to the dynamic Jahn- Teller 
effect. Using electron-phonon couplings extracted from density functional theory we calculate the 
tunneling spectrum through the C^g anionic state and find excellent agreement with measured data. 

PACS numbers: 73.63.-b, 68.37.Ef, 63.22.-m, 71.70.Ej 
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Orbital electronic degeneracy and stability of the nu- 
clear configuration are incompatible unless the atoms of 
a molecule lie on a straight line. This statement is due 
to Jahn and Teller [l| , who proved its general validity in 
1937 using group theory [2]. The theorem has important 
implications for non-linear molecules with a degenerate 
electronic ground state, particularly when the degener- 
ate electrons participate in the binding of the molecule. 
In such cases the molecule undergoes a structural dis- 
tortion which lifts the electronic degeneracy by reducing 
the symmetry of the nuclear configuration, known as the 
static Jahn- Teller (JT) effect. In certain situations sev- 
eral distortions can lower the symmetry of the molecule, 
and hence lift the electronic degeneracy. When such dis- 
torted states are also degenerate, the system will fluctu- 
ate between these equivalent configurations by quantum 
tunneling, resulting in pseudorotations 3] and restora- 
tion of the parent symmetry, denoted the dynamic JT 
effect [i,i]. 

The buckminsterfuUerene Cgo is an exceptionally sym- 
metric molecule, and its symmetry governs many of its 
physical and chemical properties. In group theory it is 
classified by the icosahedral point group Ih, the group 
with the highest number of symmetry operations in three 
dimensions [J|. The geometrical symmetry is also re- 
flected in the electronic and vibrational properties, which 
are both highly degenerate. For instance, the highest oc- 
cupied molecular orbital (HOMO) and the lowest unoc- 
cupied molecular orbital (LUMO) are orbitally fivefold 
and threefold degenerate, respectively. By adding elec- 
trons or holes to Cgo these degeneracies are explored and, 
hence, JT physics is expected to play an important role 

a. 

Indeed, the JT effect has been reported in photoemis- 
sion (PE) spectra of Ceo in the gas phase for different 
charge states of the molecule [71, Isl, |9|, llO|, lll|- Here, 
a fine structure in the intensity spectrum has been ex- 
plained as originating in the underlying electron-phonon 



(e-ph) coupled problem. Signatures of JT have also been 
found in luminescence spectra of Cgo nanocrystals [l2l |. 
A more direct consequence of the JT effect has been ob- 
served in low-temperature scanning tunneling microscopy 
(STM) images of potassium-doped fuUerenes (KaCgo) on 
Au(lll), revealing a significant static deformation of the 
molecular structure 13]. While vibronic effects have also 
been studied in the transport properties of Cgo molecules 
l4 , llSl . Ild | , the dynamic JT effect in single-molecule con- 
ductance has not been reported so far, presumably due 
to the difficulty of decoupling a molecule sufficiently from 
its environment. 

From the theoretical side there has been a substantial 
interest in describing the e-ph coupled problem in Ceo, 
primarily motivated by the discovery of the supercon- 
ducting phase of alkali doped fuUerides A3 Ceo and its re- 
lation to the "conventional" e-ph mechanism IT], [iSl, ll9| . 
An accurate calculation of the e-ph coupling constants is 
fundamental to such descriptions. In comparison with 
couplings fitted to photoemission data on C^g, it has 
been argued that couplings derived from density func- 
tional theory (DFT) are too weak 0. Contrary, DFT 
couplings were later shown to accurately account for pho- 
toemission data on C^q [10]. This apparent controversy 
calls for complementary studies of the e-ph interaction 
that put DFT-derived couplings to test. 

In this Letter, we report on the JT fingerprint in the 
conductance spectra of single Cgo molecules in a tunnel- 
ing junction. FuUerenes were decoupled from a Au(lll) 
surface by a template of organic molecules. Scanning 
tunneling spectroscopy reveals a broad sideband above 
a sharp LUMO-derived resonance. To understand this 
we calculate the tunneling spectrum by solving the JT 
problem from first principles. The vibrational modes, 
frequencies, and e-ph couplings are derived from DFT, 
and the excitation spectrum of the vibronic ground state 
is computed numerically by nonequilibrium Green's func- 
tion (NEGF) techniques as well as by exact diagonaliza- 
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FIG. 1: (color online) (a) STM image of the double rows 
of alternating TPA/Ceo on Au(lll) (/ = 16 pA, V = 580 
mV, r = 5 K). Some isolated Ceo molecules are adsorbed in 
the troughs between the rows, (b) Side view of the structural 
model along the dashed line in (a) and dl/dV spectra recorded 
on top of such a fuUerene (/set = 0.46 nA, 14ct = 2.5 V, 
T/™s=14mV). 



tion. The theoretical results explain the origin of the 
experimentally observed fine structure as due to the dy- 
namic JT effect. 

The electronic decoupling of the molecule from its 
metallic substrate is crucial for studying any JT-induced 
fine structure of the molecular resonances. To re- 
duce this interaction experimentally we use a novel 
strategy based on codeposition of Cgo and 1,3,5,7- 
tetraphenyladamantane (TPA) on Au(lll), which results 
in spontaneous formation of nanostructures where the 
Ceo cages are lifted away from the surface by support 
from TPA molecules 22|, |23| . One of the observed struc- 
tural motifs are double rows of alternating TPA and Cgo, 
cf. Fig. [IJa). This dielectric template yields troughs for 
single Cgo molecules to be trapped between neighbor- 
ing double rows as sketched in Fig. [IJb). Differential 
conductance {dl/dV — V) spectra taken over the center 
of these molecules reveal strong non-linearities around 
-2.2 V, 1.3 V and 2.5 V, associated with the HOMO, 
LUMO and LUMO-fl derived resonances, respectively 
(Fig. mb)). The corresponding large gap of ^ 3.5 eV 
is a clear sign of a weakly interacting molecule with its 
surrounding [23|. Even more striking is the observation 
of a very sharp LUMO-derived resonance (FWHM r^ 60 
meV), reflecting a long lifetime in the order of 10 fs of the 
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TABLE I: Partial electron-phonon coupling constants 
X,./N{0) = 2/9 J2i,j \Mijf/hLO^ of the Ag and Hg intramolec- 
ular modes for Ceo', cf. Eq. (10) of Ref. [H, with iV(0) being 
the density of states. The sum of all modes gives a coupling 
strength of X°^^/N{0) = 76.1 meV S in good agreement 
with previous calculations [ISl . Il9l . l20l . l2ll | . 



tunneling electron. Hence, this system is ideal to study 
e-p h couplings. However, instead of single vibronic peaks 
[l6|, we observe a broad sideband at ~ 230 meV above 
the LUMO energy, which lies well outside of the 200- 
meV-wide vibrational spectrum of Ceo. 

Tunneling electrons through an isolated Cgo molecule 
will probe the excitation spectrum as resulting from a 
dynamic JT interaction with the vibrations. Near the 
LUMO-derived resonances in focus here, we exclude the 
HOMO and LUMO-l-1 states from our treatment since 
these are well separated in energy. Further, group theory 
gives that only the Ag and Hg intramolecular phonons 
couple to the ii„ LUMO states 6]. Thus, the JT problem 
for CgQ is described by the following Hamiltonian 

3 42 



H = eo^ c\ci + ^ hwublK 
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+Y.Y.Y.^'^lA^M+b<^)^ (1) 

IV=1 2=1 J = l 

where c\ (q) is the one-electron creation (annihilation) 
operator corresponding to one of the three degenerate 
single-particle LUMO states \i) and hi, {b^) is the bosonic 
creation (annihilation) operator of each of the 42 (8 five- 
fold degenerate Hg and 2 non-degenerate Ag) vibrational 
modes v. The parameters Eq, uj^, and M^, represent the 
bare LUMO energy, the vibrational frequency of mode v^ 
and the coupling constant for the scattering of an electron 
from state \i) to state \j) under creation or annihilation of 
a phonon in mode v^ respectively. These parameters (see 
Tab. HI) are obtained from DFT calculations 2A IH, Ilol 
applying the scheme described in Ref. [27] to the neutral 
,21.] . isolated Ceo molecule. 

In our STM configuration the molecule under investi- 
gation is much less coupled to the tip than to the sub- 
strate, i.e.. Ft <S Fs in terms of the tip Ft and substrate 
Fg tunneling rates. Under this condition the electron oc- 
cupancy of the molecular states is effectively in equilib- 
rium with the substrate and the differential conductance 
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FIG. 2: (color online). (a) Theoretically computed JT- 
induced fine structure in the LDOS of the LUMO via the 
coupling to Ag and Hg modes, shown for two different in- 
trinsic level broadenings: F = 10 meV (black line) and F = 
60 meV (thick blue line), (b) Analysis of the position of the 
phonon sidebands including only one mode at the same time, 
calculated with both the SCBA (black lines) and LS (dashed 
red lines) methods (F = 10 meV). For clarity the spectra are 
off-set and shifted to place the main peak at zero energy. The 
vertical dotted lines indicate the fundamental frequencies of 
the Ag and Hg modes. 



at a bias voltage Vs is to a first approximation propor- 
tional to p(/is + eVs)Tt, where p is the local density of 
states (LDOS). This takes into account both coupUngs 
to the leads as well as to the vibrations. We have con- 
sidered two different methods for calculating p, namely 
exact diagonalization with the Lanczos scheme (LS) as 
well as NEGF within the self-consistent Born approxi- 
mation (SCBA) m. 

Figure ^a.) shows the computed JT spectrum from 
SCBA for two different values of the total electronic 
broadening F = Ft + F^ «:! F^ (FWHM). The spectrum 
is composed of a main peak followed by a series of side- 
bands at higher excitation energies. Depending on the 
broadening, the LUMO spectrum might appear as one 
main peak with a single sideband separated by 230 meV. 
A similar structure has also been reported in photoemis- 
sion [71, la y, lldl ■ It is important to realize that the fine 
structure in the LDOS does not relate in a simple way to 
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FIG. 3: (color online). Direct comparison between experi- 
ment (red circles) and theory (black line) for three different 
molecules: (a) /set = 0.3 nA, Vset ~ 2.5 V, Kms ~ 7 mV, 
Fs = 60 meV. (b) /set = 0.52 nA, Ket = 2.5 V, Kms = 3.5 
mV, Fs = 60 meV. (c) /set = 1.2 nA, V^et = 2.5 V, Kms = 7 
mV, Fs = 50 meV. Panel (a) includes calculations with scaled 
e-ph couphng constants (dashed blue lines) corresponding to 
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the fundamental frequencies of the Ag and Hg modes (in 
the range 32-195 meV as indicated with vertical dotted 
lines in Fig. ^. Hence, for JT systems it is not possi- 
ble to directly relate vibronic structure in spectroscopy 
with molecular frequencies. Another important fact is 
that the vibronic spectrum is threefold degenerate, i.e., 
the parent symmetry of the LUMO is restored by the Ag 
and Hg modes as expected from the dynamic JT effect 
[29[. Hence, the molecule is not statically distorted. 

To gain an understanding of the contributions from 
the different vibrational modes to the total spectrum, we 
show in Fig. [2{b) the LDOS as resulting from a calcu- 
lation with the coupling to only one type of mode at a 
time. Although the sideband structure is much simpler, 
one observes that the peak separation for all Hg modes is 
larger than the phonon energy. This nontrivial behavior 
is generic to the Tiu®hg JT problem [6]. Figure[2l^b) also 
shows that, except for the more strongly coupled Ag{2) 
mode, SCBA and LS yield essentially the same results 
for the LDOS due to the weak e-ph couplings. 

We next turn to a comparison with the experiment. 



The measured tunneling spectra on different molecules 
exhibit certain variations in peak position and structure 
around the LUMO resonance. This variation is likely 
due to slightly different environments, and possibly also 
molecular orientation, of the selected Ceo species. In 
Fig. [3] we compare the dI/dV-specti& of three different 
molecules with theoretical spectra by applying an appro- 
priate substrate coupling F. The relative peak height as 
well as the position of the sideband with respect to the 
main peak are in very good agreement. This fact sup- 
ports the interpretation that these spectra indeed exhibit 
the free-molecule JT effect. 

Alternative mechanisms, such as the weak spin-orbit 
coupling in carbon [J] or the Stark effect from the applied 
voltage, cannot account for the 230 meV splitting seen 
in Fig. [3] Although details in the weak electronic cou- 
pling of the Ceo molecule to its environment may have an 
influence in the experimentally recorded spectra, we be- 
lieve that this effect solely cannot produce such an agree- 
ment as the JT theory presented above for a decoupled 
molecule. Furthermore, from the different timescales 
of electron tunneling events and residence times in the 
molecule, we can also exclude heating effects. Based on 
the NEGF-SCBA calculations we find that heating would 
be important if the current through the LUMO states is 
raised by two orders of magnitude. 

In Fig. E^a) we also show how the calculated spectrum 
changes by a scaling of the DFT-derived e-ph couplings 
corresponding to X^ -^ 1/2 A^°^ and X^ -^ 2 A^*-^^. As 
seen, these scalings decimate the agreement with experi- 
ment in terms of peak ratio and separation, thus pointing 
towards reasonable values of the e-ph coupling strength 
within DFT. 

In summary, we have explained the dynamic JT ef- 
fect as it is revealed in electron tunneling spectroscopy 
through the Cgg anionic state of sufficiently isolated 
molecules. The experimental realization of this situation 
was based on a novel experimental preparation procedure 
of codeposition of Ceo and TPA molecules on Au(lll), 
where Ceo molecules on top of the double row template 
are exceptionally decoupled, geometrically and electron- 
ically, from the metallic substrate. Low-temperature 
STM dl/dV spectra recorded under such conditions dis- 
play a delicate fine structure around the LUMO-derived 
resonance, that closely resembles the theoretically com- 
puted JT spectrum. This quantitative agreement further 
supports that the calculated e-ph couplings from DFT 
provide an accurate description for the vibrational inter- 
actions in the molecule. 
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